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Abstract 

; We prove that there exists a strong solution to the Dirichlet boundary value problem 

for the steady Navier-Stokes equations of a compressible heat-conductive fluid with large 
external forces in a bounded domain f2 C M d (d — 2, 3), provided that the Mach number 
is appropriately small. At the same time, the low Mach number limit is rigorously verified. 
The basic idea in the proof is to split the equations into two parts, one of which is similar 
to the steady incompressible Navier-Stokes equations with large forces, while another part 
corresponds to the steady compressible heat-conductive Navier-Stokes equations with small 
forces. The existence is then established by dealing with these two parts separately, estab- 
lishing uniform in the Mach number a priori estimates and exploiting the known results on 
, the steady incompressible Navier-Stokes equations. 
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1 Introduction 



This paper is mainly concerned with the existence of strong solutions to the steady Navier- 
Stokes equations of a compressible heat-conductive fluid in a bounded domain Q C M rf (d = 2, 3) 
with large external forces: 

div(^u) = 0, 

< gu ■ Vu + Vp = divS(Vu) + g{ + g, (1.1) 
^^U' V9 +pdivu = kA0 + \&. 

Here g denotes the density, uGl li the velocity, the temperature, p = RgQ the pressure with 
R > being the gas constant, c v > is the heat capacity at constant volume; f is the density 
of external body force and g is a given external force. The stress tensor § and the dissipation 
function ^ are defined by 

S = 2fiD(u) + Adivul, * = 2fiD(u) : D(u) + A(divu) 2 > 0, 



* Corresponding author. E-mail address: douchangsheng@163.com (C. Dou), jiang@iapcm.ac.cn (S. Jiang), 
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D(u) = (Vu+Vu')/2 is the deformation tensor. The viscosity coefficients fi, A satisfy 2/j>+d\ > 
and fi > 0, k > is the heat conductivity coefficient. Moreover, the total mass is prescribed: 

gdx = M > 0. 

We impose that the velocity u satisfies no-slip boundary condition and the temperature 
is a constant $0 on the boundary of , i.e., 

u = 0, e = i? on an. (1.2) 

In the last decades, the steady compressible heat-conductive Navier-Stokes equations have 
been studied by many mathematicians and there are a lot of results on the existence in the 
literature, here we recall some of them for both small and large external forces which are related 
to our study in this paper, and we refer to the monograph |20| for more details. When external 
forces are sufficiently small, Matsumura and Nishida in 1982/83 proved the existence of a solution 
with potential forces near a rest state [141 H5j . while Valli and Zajackowski |23|. [25] used the 
existence of global non-stationary solutions to get the existence of stationary solutions. Later, 
Valli [24J showed the existence of stationary solutions in the general case by using an idea of 
Padula [21] to decompose the equations into two parts that are governed by the Stokes equations 
and a transport equation, respectively. Beirao da Veiga |T] obtained more general existence 
results in the L p -setting by decomposing the equations into three parts that are governed by the 
Stokes equations, a transport equation and Laplace's equation, respectively. In 1989, Farwig 
[5] showed the existence of solutions to the steady compressible heat-conductive Navier-Stokes 
equations for small forces with slip boundary condition. 

When external forces are of arbitrary size, the existence of strong solutions was proved 
in [T5J [16] for the case of potential forces. When the equations of state and the viscosity 
coefficients satisfy certain (growth) conditions, Novotny and Pokorny |17[ [T8] showed that weak 
or strong solutions to the steady compressible heat-conductive Navier-Stokes equations exist. 
Unfortunately, their results exclude the case of ideal polytropic gases, for which the existence of 
strong solutions, to our best knowledge, still remains open. 

The aim of the present paper is to establish the existence of strong solutions to the steady 
compressible heat-conductive Navier-Stokes system (II. ip without any smallness assumption on 
the external forces f and g, when the Mach number is small. 

We mention that in the isentropic flow case, the existence of weak solutions or strong solutions 
under small Mach number for large external forces has been extensively investigated. Lions [13] 
first proved the existence of weak solutions under the assumption that the specific heat ratio 
7 > 1 in two dimensions and 7 > 5/3 in three dimensions. The restriction on 7 actually comes 
from the integr ability of the density g in L p , and in fact, the higher integr ability of g has, 
the smaller 7 can be allowed. In [20] Novotny and Straskraba showed the existence of weak 
solutions for any 7 > 3/2 if f is potential and g = 0. By deriving a new weighted estimate of the 
pressure, Frehse, Goj and Steinhauer [6], Plotnikov and Sokolowski |22] established an improved 
integrability for the density under the assumption of the L 1 -boundedness of gu 2 which was not 
shown to hold unfortunately. In 2008, Bfezina and Novotny [3] was able to prove the existence of 
weak solution to the spatially periodic problem for any 7 > (3 + \/41)/8 when f is potential and 
g = 0, or for any 7 > (l + \/l3)/3 ~ 1.53 when f, g 6 L°°, without assuming the L 1 -boundedness 
of gu 2 , by combining the L°°-estimate of A _1 P with the (usual) energy and density bounds. 
Then, in the framework of [3], Frehse, Steinhauer and Weigant [3 |8] established the existence 
of weak solutions to the Dirichlet boundary value problem for any 7 > 4/3 in three dimensions 
and to the spatially periodic or mixed boundary value problem for 7 = 1 (isothermal flow) in 
two dimensions. Recently, Jiang and Zhou \10\ fTT] proved the existence of weak solutions to the 
spatially periodic or Dirichlet boundary value problem in M 3 for any 7 > 1. The existence of 
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strong solutions was shown by Choe and Jin [3] when the Mach number is small, by exploiting 
the known results for the incompressible steady Navier-Stokes equations. 

Now, we rewrite (jl.ip in the form of the Mach number. After scaling and a straightforward 
calculation we obtain the following dimensionless form of the steady full compressible Navier- 
Stokes equations: 

div(£>u) = 0, 



gu.Vu+4 = divS(Vu) + gf + g, 
gu ■ VG + p divu = nA@ + e 2 ^, 

where e is the Mach number. 

Since the total mass of the fluid is given, we impose the condition 



;i.3) 



g := — / g{x)dx > 0, 



which can be renormalized to g = 1 without loss of generality. Similarly, we also assume that 



9 = 1, R 



1, #0 = 1. 



To show the existence, we take the transformation 

g=l + ep, = 1 + e6 
to rewrite the system (jl.3p in the form: 
divu + ediv(pu) = 0, 

(1 + ep)(u ■ Vu) + £±IpVP + = divS ( V u) + (1 + ep)f + g, 

e(l + ep)u • V6 + divu + (ep + ed + e 2 p9)divu = eK,A6 + e 2 ^, 
with boundary conditions 

u = 0, 6 = ondn. 
Now, we state the main result of this paper. 



:i-4) 



;i.5) 



;i.6) 



Theorem 1.1. Let f , g G H 2 (9). Then there is an eo depending on ||(f, g)||ij 2 an d ^> such 
that for any e £ (0, eo), there exists a solution (p e , u e ,9 e ) G H 2 x (H 3 n Hq) x (H 3 n i?o) *o f/ie 
boundary value problem (]1.5I) . (jl.6p . satisfying 



lim inf ||u e 

e^0(7,PeL 



^113 + 11^112 + 



P £ + 



P 



where (U,P) G L := {(Z7, P) G (H a C\Hq) x PT 3 | (U,P) is a solution of the incompressible steady 
Navier-Stokes equations (|1.7p with external force f + g}, i.e., 



U-VU- pAU + VP = f + g, 
div U = 0, 



u = o on an, 



1.7) 



Pdx = 0. 



Remark 1.1. If considering the existence of spatially periodic solutions to (|1.5p in a periodic 
domain, we can also obtain a existence result similar to Theorem \l.l\ 



The system (|1 . 5[) is complicated mixed-type nonlinear equations containing such structures as 
elliptic and hyperbolic systems, for which the usual approach of the fixed point arguments used 
to prove the existence of classical solutions requires the smallness of data. To show Theorem 
1.11 (the existence for large data), we split the system (jl.5p into two parts, one of which is 



similar to the steady incompressible Navier-Stokes equations with large force f +g, while another 
part corresponds to the steady compressible heat-conductive Navier-Stokes equations with small 
force ef, provided the Mach number e is small. Then, as noted in [3], we modify and combine 
elaborately the arguments in |9j where the existence of strong solutions to the incompressible 
Navier-Stokes equations for large forces was presented, and in [5] where strong solutions of the 
compressible viscous heat-conductive equations with small forces were dealt with, to establish 
Theorem 11.11 

Compared with the isentropic case studied in [3], due to presence of the energy equation, the 
main difficulties here lie in the existence of weak solutions to the approximate linearized system, 
dealing with the coupling terms between the velocity, density and temperature, and deriving 
the uniform estimates in a bounded domain, for example, how to control the energy norm 
|| u — U\\a + H77H2 + ||#||3 uniformly in e under the no-slip boundary condition. To circumvent such 
difficulties, we take the transform of g = 1 + ep, O = 1 + e8 for the system (|1.3j) . instead of the 
transform (g = 1 + e 2 p, = 1 + e 2 9) used in [3], and utilize the lower order terms to control the 
higher order terms. More precisely, the main steps of the proof are the following: First, we apply 
Lax-Milgram's theorem to get the existence of weak solutions to the regularized elliptic equations 
f)2.14[) ([2.16|) of the transformed linearized equations. Then, we exploit the uniform estimates 
and a compactness argument to get the existence of a weak solution to the transformed linearized 
system (j2. 12|) . We can easily verify that the weak solution to the transformed linearized system 
(12.12p is indeed a weak solution to the linearized system (I2.10p . And this fact together with 
the uniqueness of weak solutions to the system (|2.10p gives the existence of a weak solution to 
the approximate linearized equations (12.10P in Section 2.2. Second, we exploit the property of 
the momentum equations and the regularity of the Stokes problem to establish the estimates 
of \rj + 9\, which, combined with an estimate for 8, implies the boundedness for rj. Due to 
presence of boundary here, some difficulties involved with controlling the boundary terms arise. 
To overcome such difficulties, the crucial step is to get a .£f 2 -bound of divu near the boundary, 
for which we shall adopt the local isothermal coordinates used in |231 125] . This strategy has 
also been used in [2j [12] to study the low Mach limit of the compressible Navier-Stokes with 
non-slip boundary condition. Then, summing up all the estimates for (v, 9) and 77, we can 
establish the desired a priori uniform in e estimates in view of the smallness of e (see Section 3). 
Finally, we apply the Tikhonov fixed point theorem to obtain the existence of a strong solution. 
Moreover, with the help of the uniform a priori estimates, one can take to the limit to show 
the incompressible limit. We point out that due to the splitting, we have to impose that the 
energy equation should not possess a heat source. Furthermore, it is still open whether a strong 
solution of the steady incompressible or compressible Navier-Stokes equations is unique. 

The paper is organized as follows. In the next section, we will prove the existence of weak 
solutions and regularity to the linearized incompressible and compressible problems. Section 3 is 
devoted to establishing the existence for the nonlinear problem. Finally, the incompressible limit 
of solutions to the steady compressible heat-conductive Navier-Stokes equations is presented in 
Section 4. 

Notations: We denote by 1? the Lebesgue space L 2 (VL) with norm || • ||o, by H m the Sobolev 
spaces H m (£l) with norm || • ||m* Define the spaces 

H m = {p€H m \ j^p(x)dx = 0} 1 F 1 >(J = {uG J ff 1 |divu = o}, fl£. := H m n H^. 

We denote by H^ 1 the dual space of H$ with the dual product (•, •) and the norm || • = 
sup|| fe || 1=1 \(-,h)\. We shall use the abbreviation: 

J ■ dx := J ■ dx. 
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2 Existence of solutions to the linearized problem 



We first split the system (jl.5p into two parts, so that one part looks like the incompress- 
ible Navier-Stokes equations, while the other part behaves like the compressible Navier-Stokes 
equations. More precisely, let (U, P) and (v , 77) be the solutions to the following systems, respec- 
tively: 



U ■ W + v • W - fiAU + VP = f + g, 
div U = 0, 



P = 0, 



(2.1) 



U = on dO, and 



P = 0; 



and 



U • Vr/ + 



divv 



-v • V?? — j^divv — ediv(P([7 + v)), 



U • Vv - /jAv - + A)Vdivv + V? ? + V6> = eF - v • Vv - (9Vr/ - ??V(9, 

i . e (2.2) 

U-V9- kA9 + = eG - v • V0 - r/divv - fldivv, 



0, 



9 = on <9S1 and 



7? = 0, 



where the new force F and heat source G are defined by 

F = (eP + r?)f - (eP + r?)(C7 + v) • V{U + v) - #VP - PV9, 
G = ^ - (eP + rj)(U + v) • V9 - (eP + ry^divv - Pdivv. 

It is clear to observe that u := U + v, p := eP + r] and 9 are a solution to (jl.5p . Thus, we can 
obtain a solution of the system (jl.5jl if we can solve the systems (|2.ip and (12. 2h . First, we will 
give the existence of weak solutions to the linearized incompressible problem (|2,ip and derive 
a priori estimates of higher order derivatives of the unknowns (U,P). Then, we shall show the 
existence of weak solutions to the linearized compressible problem (|2.2p and establish uniform 
estimates of higher order derivatives of the unknowns (77, v, 9). 

In what follows, we assume that meas(f2) = 1 without loss of generality. 



2.1 Linearized incompressible equations 

Let U and v be given functions satisfying U G H 4 n H^^ and v G H s n Hq. At first, we 
consider the linearized equations to (|2.ip for given U and v as follows. 



(U + v) • W - fiAU + VP 
divt/ = 0, 



h. 



P = 0, 



(2.3) 



U = on <9S1 and 



P = 0. 



where h = f + g. 

The problem (|2.3p is a Stokes problem which is solvable for arbitrarily large forces. In fact, 
(|2.3p can be solved by using the Lax-Milgram theorem for small v, and we can obtain the 
following existence result, the proof of which can be found, for example, in 0], and is therefore 
omitted here. 

Lemma 2.1. Let h G H^ 1 . There exists a constant ao depending on fi and Q,, such that if 
1 1 v 1 1 3 < ao, then there exists a weak solution (U,P) G Hq x H° of (j2.3j) . satisfying 



\U\\x < Co\\h\\-i, 
5 



(2.4) 



iiPiio^diihiua + Hhiu), (2.5) 

where Cq and C\ are positive constants which depend only on O, \x and a$. 
As for the regularity of solutions, we consider the Stokes equations: 

—fj,AU + VP = h - (U + v) • VP, 
divP = 0. 

Then we can derive the following estimates by employing bootstrap arguments similar to those 
in @]. 

Lemma 2.2. Let h G H m , U G P m+1 n P 1 , m = 0,1,2, and v be the same as in Lemma 
\2.1[ There are positive constants C 2 , C3 and C4, depending only on VL, [i and ao, such that if 
U G Hq satisfies the inequality (|2.4p , i/ien 

||^||2 + ||VP||o<C 2 ||h||o(||h|| + l) 4 . (2.6) 



LfUeH 2 nH^ satisfies ([M]), i/ien 

||C/|| 3 + ||VP||i ^CgllhUidlhUi + l) 8 , (2.7) 
and i/P G P 3 n P 1 satisfies (f2"77|) . £/ien 

II^IU + ||VP|| 2 < C 4 ||h|| 2 (||h|| 2 + l) 12 . (2.8) 
Let f,g £ P 2 (f2), then it is obvious that h G P 2 (f2). We define a function space Kq by 

P := {Pe<P): ||P||i < Cl||h||i, \\U\\ 2 < C 2 ||h|| (||h|| + l) 4 , 

\\U\\ 3 < CallhUidlhll! + l) 8 , \\U\\ A < C 4 ||h|| 2 (||h|| 2 + l) 12 }. (2.9) 

Thus, by Lemma 12.21 we see that the solution P of the system fj2.3j) also lies in Kq for any 
given U G Kq, since the constants C\, ■ ■ ■ , C4 do not depend on U. 

2.2 Linearized compressible equations 

Let (P, v, 9) G (P 4 n Hl a ) x (P 3 n P 1 ) x (P 3 n P 1 ) be given functions, and (P, P) be the 
solution of (|2.3p established in Section 2.1. Next, we give the existence of weak solutions and 
derive some a priori estimates for solutions to the linearized equations of the system (12. 2D . For 
simplicity, we only consider the three-dimensional case. As aforementioned, we shall apply the 
Tikhonov fixed point theorem to show the existence of steady strong solutions to (I1.5P . To this 
end, for given (v, 9) G (P 3 fl Pq) x (P 3 nPo), let (??, v, 9) be the unique solution of the following 
linearized system of (12. 2p the existence of which will be shown below: 

' divv 
U ■ Vt] H h v • Vr/ + r/divv = -ediv(P(P + v)), 

U • Vv - /iAv - CVdivv + V?? + — + 6Vt/ + r/Vfl = eP - v • Vv, (2.10) 

P • V0 - nA9 + + r/divv = e(5 - v • V0 - fldivv, 
e 

where the new force F and heat source G are defined by 

F = (eP + r?)f - (eP + rj)(U + v) • V(P + v) - 9VP - PV9, 
G = § - (eP + rj)(U + v) • V0 - (eP + r/)0divv - Pdivv, 
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with § = 2fiD(U + v) : D(U + v) + A(div(£7 + v)) 2 and C = fi + A. 
Thus, for given U, v and 0, we can construct a map TV: 

N{U,v,0) := (I7,v,0). 

And, we have to show that N maps some space into itself and is weak continuous to get a fixed 
point of the mapping N. 

In order to obtain the existence of weak solution of (|2.10j) . we set 

F' = ePi - eP(U + v) • X7(U + v) - 6VP - PX78, 
& = § - eP{U + v) • V0 - ePOdivv - P divv. 

So, F = F' + nf + rj(U + v) • V(C/ + v), G = & - 7](U + v) • V0 - r^divv. 
2.2.1 Existence of weak solutions 

Lemma 2.3. Let F',G' E H , f € -ff 2 , and (U,P) be a solution of (|2.3p established in Lemma 
\2.1\ If 1 1 v 1 1 3 + ||0||3 is sufficiently small, then there exists a unique weak solution (n,v,0) £ 
H° x Hq x Hq to the equations f)2.10|) boundary condition (|2.2j) 1 . 

Proof. First, the momentum equations (|2.1U|) 2 can be rewritten as 

T) + + erfi + eA _1 div(C7 • Vv) - e(/u + C)divv + e 2 A _1 div n(U + v) • V(U + v) 

- e 2 A~ 1 div(r ? f) = eA _1 div(eF' - v • Vv). 

We add ()2.11j) to ()2.10p i and ()2.10j) ^. respectively, and rewrite the system (|2.10p as the following 
equations: 

' divv ~ i 

U • V?? H h v • Vr/ + r/divv + rj + 6 + er]9 + eA _i div(C7 • Vv) 

e 

-e(fi + C)divv + e 2 A~ 1 div rj(U + v) • V(U + v) - e 2 A- 1 div(??f) 
= -ediv(P(C7 + v)) + eA^div^F' - v • Vv), 

(/•Vv- /iAv - CVdivv + V?? + V61 + OVrj + ryVfl - er/f + er](U + v) • V(Z7 + v) 

e 

= e i?' _ v • Vv, 

divv ~ i 

U -V6 - kA6 H h r/divv + er/(£/ + v) • V0 - e??0divv + + + er/0 + eA _1 div(C/ • Vv) 

e 

-e(/i + C)divv + e 2 A _1 div rj(U + v) • V(U + v) - t^A^divfaf) 
= eG' - v • V0 - 0divv + eA^div^F' - v • Vv), 



(2.11) 



with boundary conditions 

v = e = o on on. 

Consider the regularized elliptic equations of (12.12P as follows: 



(2.12) 
(2.13) 



-5Arj 5 + U-Vn s + 



divv 



+ v • Vr/ 5 + 7/divv + n s + 9 s + erf~B + eA _1 div(i7 • Vv 5 



e(n + Odivv 5 + e 2 A _1 div \n 5 {U + v) • V(U + v)l - e 2 A^div^f ) 



-ediv(P(C7 + v)) + eA -1 div(eF'" - v" 5 • Vv 



(2.14) 



U ■ Vv 5 - fiAv s - CVdivv 5 + Vr ^ + W * 



+erf(U + v) • V(U + v) = eF' S - v 5 • Vv, 



(2.15) 



u ■ ve d - kA6 6 + 



divv 5 



+ r/divv + erf(U + v) • V0 - e^fldivv + rf + 6> 5 + erf 



+eA~ 1 div([/ • Vv 5 ) - e(/u + C)divv 5 + e 2 A- 1 div \r] S (U + v) • V(U + v) 
-e 2 A~ 1 div(r/' 5 f) = eG' - v" 5 ■ V6 - 5 divv + eA -1 div(eF' - v • Vv), (2.16) 
with boundary conditions 



^ = 0, vW 
an 



on dQ. 



(2.17) 



Here n is the outer normal vector. 

The system (|2.14p - (|2.17p in variational form reads as: Find (r] 5 ,v 5 ,9 5 ) £ H 1 x Hi x H^, 
such that 



£(77*y,0*;77,v,0) 



[-ediv(P(£7 + v)) + eA" 1 div(eF / - v • Vv)]r/ + (eF' - v • Vv)v 



+ [eG' - v • V0 - 5 divv + eA _1 div(e.F' - v • Vv)]9dx, 



(2.18) 



where 



B{rj S ,v S ,8 5 ;ri,v,9) ■= 5 J V? / • V !Z dx - yV ■ V !?) r ?' 5 + • V X) • v 5 + ([/ • V9)9 s dx 

- J (v • Vr^ifdx + y (?7 5 + 6» 5 )(r? + 0) + /xVv* 5 : Vv + C(divv <5 )(divv) 
+/tV6» 5 • VMc + /" jer/^ + eA -1 div(?7 • Vv 5 ) - e(/i + C)divv' 5 
-e 2 A- 1 div(r/' 5 f) + e 2 A _1 div ^([7 + v) • V(U + v) }(?? + 0)dx 
er7 5 (f7 + v) • V(J7 + v) - er] s {\ ■ v - rfO&vvdx 



er] S (U + v) ■ V0 - er/^divv + r/divv 
+- y (divv 5 )?? - divv(r/ + 5 ) + (divv 5 )0cix, 



(2.19) 



for any (rj,v,9) e H 1 x Hq x Hq. 

We shall apply Lax-Milgram's theorem to show the solvability of the variational problem 
(|2.18p . We first show the coerciveness of the bilinear form in (j2. 19|) . To this end, we take r] S = 
V, V 5 = v, 9 s = 9 in (I2.19p . integrate by parts and use Poincare's inequality ||0 a ||o < ^/co||V0 ||o 
to see that 



B(rf ,v ,9°;r] 

> y (5\Vri s \ 2 + (rf + 9 5 ) 2 + /x|Vv 5 | 2 + ((divv 5 ) 2 + K \V9 5 \ 2 )dx - J 



\rj 5 \ 2 )dx 



< I u li|Vv 5 || 2 - ||0|| 2 ||r/|| 2 - ||0|| 2 ||divv 5 || 2 - ||divv|| 2 (||^|| 2 + ||0 5 || 2 ) 



-e(||f|| 2 + ||C/ + v||^ + ||[/ + v 



13 + l|v|| 3 | 



-e(||f || 2 + \\U + nl)\\v s \\ 2 - e(\\U + v|| 2 ||0|| 3 + ||v|| 8 ||0|| 2 ) 



2(111^112 + wnt) + nio + raiivvig + Oi + crnio + iiv^h,-; 

-r + vi||^||g + ||^||2 + e ||f|| 2 (||^||2+||^||2)] 



a<5 m2> 



3<5i|2 



<5i|2 



?'llg 

\2||/i<5||2 



J\\2 



> 



2e\\e\\ 2 - e(v + O 2 ~ \\Oh 
% + llvllsl 



2co + K 

-2||?7 + v||2 + ||;7 + v 
2 



M -2e||C/||^-ec (||f|| 2 + ||C/ + v' 5 ||| 



L2 



C o||divv|| 2 -e Co (||C/ + v 5 || 2 | 



||divv|| 2 - (e + e 2 V^)(2||f|| 2 

n s \\ 2 + s\\v v 5 \\ 2 

l|Vv 5 ||2 + ( C - ~ 

3 + ||v|| 3 W2 + || 



2 -e)||divv 5 [|g 



+i + (^ + C) 2 + ||f|| 2 
> S\\Vr, s \\ 2 + C'(\\ V 5 \\ 2 



||V6l'||g 



+ v 



<5i|2 



+ \\o s \\i), 



(2.20) 



where the constant C depends only on $7, fi, A, k, Co, ||0||2 and 1 1 v 1 1 3 . Notice that here the 
smallness of e, ||#|| 2 and 1 1 v 1 1 3 are necessary. In fact, we have used the following estimate in the 
second inequality of (|2.20|) : 



f n v s + 9 5 ) 2 + -\v9 8 \ 2 )dx> [ Un 5 + e 5 f + —\e & \ 2 )dx 

J 2 J 2co 



2co ,/? + 



2cq + K 



2cq + K 



\ v s \ 2 y x 



> 



2cq + K 



hi 



Obviously, the bilinear form in (|2.19|) is continuous on H 1 x Hq x Hq. Thus, there is a 
unique solution , v s ,0 s ) of (|2.14p — ()2.17j) satisfying the variation form (|2.18p for all (77, v, 9) E 
H l x Hq x Hq. Moreover, taking e so small that e < y, we obtain 



*liv^||g + c(||i^||g 



5h2 



,<5||2 



\\e 5 \\l" 



< ||ediv(P(C7 + v)) + eA^ 1 div(eF' - v • Vv)||g + \\eF> - v • Vv||g 
+ \\eG' - v • V0 - # 5 divv + eA~ 1 div(eF' - v • Vv)||g 

for some positive constant C. 

By Rellich's compactness theorem, there is a subsequence rf™ G H°(Tl), v*' 1 G Hq(TI), 
9 s " G H^n), and r] G £f°(Q), v G flj(n), G flJ(J2) satisfying that 

rf n converges weakly to rj in H°(£l) 



and 



v «5 n ^ qSu converges weakly to v, 9 in Hq (Q) 
v Sn , 9 Sn converges strongly to v, 9 in H°(Q), 



respectively, as n — > 00. And it can be easily validated that (t],v,9) is indeed the unique weak 
solution of (12.12j) . Due to the lower semicontinuity of the iJ m -norm (to = 0, 1) and the estimate 



o5 n 112-1 



mv dn u+(\\v dn u+\\^ n \\t+w uj 

< ||ediv(P(£7 + v)) + eA^ 1 div(eF' - v • Vv)||g + \\eF' - v • Vv||g 
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+ \\eG> - v • V0 - 0divv + eA -1 div(eF' - v • Vv)[|g, 



we find that 



MI0 + H1 + 



<||ediv(P(^ + v)) + eA" 1 div(eF' - v • Vv)||g + \\eF> - v • Vv||g 



+ \\eG' - v • V9 - 6>divv + eA _1 div(eF' - v- Vv 
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On the other hand, it is easy to verify that a weak solution of (|2. 12[) is also a weak 
solution of (I2.10p . Moreover, a weak solution of (|2,10p is unique. In fact, assuming that 
(?7i,vi,#i), (?72,V2,#2) £ H° x Hq x Hq are two weak solutions to the problem (|2.10p with 
boundary condition (j2.2|) 1? and letting fj = rji — 772, v = Vi — v 2 , 6 = #1 — 8 2 , we find that 
(77, v, 9) is a weak solution of the following boundary value problem: 



U ■ Vfj + 



divv 



+ v • Vr/ + r/divv = 0, 



V?7 + V9 

U ■ Vv — /iAv — (Vdivv H h OVrj+fjVe = e(rjf + fj(U + v) • V(J7 + v)), 

- divv £ - - ( 2 - 21 ) 

U -V6 - kA6 H h fjdivv = e(-fj(U + v) • V9 - rjOdivv), 

v = 0, 9 = on ffi and / rj = 0. 



Then, we test the equations fl2.21.P p (|2.2ip 2 . (|2.2ip 3 with 7],v,6 respectively to deduce that 

/i|Vv| 2 + C|Vv| 2 cte+ / K\V9\ 2 dx (2.22) 



77V 



f + (U + v) • V(C/ + v) - J]6 (U + v) • V0 + 0divv 



— divv I 77 1 2 + 770divv0 + r?0divvdx 



l^lb + llvl 



+ l|v|| 3 (l + 



< (\\r}\\o + \W\\t 

+^(||?7|| 2 + ||v|| 2 )(||f||2 + ||C/||2 + ||v||2). 

And according to the inhomogeneous Stokes equations: 



-/iAv H 

e 

divv = divv, 
v = 0, 



-E7 • Vv + CVdivv + e 



77T + r7(C/ + v) • V(U + v) - (9Vrj + J]Ve), 



we have 



e ll v lli + II 7 ? + 0||o 
< eC{||C/|| 2 ||v||a + ||v||i + e||r7|| 



If|l2 + (ll^ll2+l|v| 



+ IMIol 



(2.23) 



By combining (|2.22j) and (|2.23p . making use of Poincare's inequality and the smallness of e, 
|v|U and II ^ II3, we see that 



Mlo + ||v||i + ||0||i <0, 



(2.24) 



which implies rj\ = r/2, vi = v 2 , 0\ = 62- Hence, there is a unique weak solution (77, v, 0) G 
H° x Hi x Hi of the problem ([2l0]) . □ 
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In order to get higher order uniform in e estimates of the (rj,v,9), we give another bound 
about ||v||i and ||0||i. 

Lemma 2.4. Let (U,P) be the solution of (|2.3p given in Lemma{2Jl Let F,G £ H~ l . Then 
we have the following uniform in e estimate: 



h\\\ + \\e\\\ <C 5 [||v|| 3 |W|g + e 2 (||^||-i + \\G\tJ + e||P|| 2 (||tf|| 2 + ||v|| 2 )||t/|| 

+ \\nf+u\i\\m+\\ntm\i + NI?) 

where the constant C5 > is independent of e. 



(2.25) 



Proof. Multiplying (|2. 10|1 -[_ . (|2.10p 2 and ()2.10P o by 77, v and 9 in L 2 respectively, and summing 
up the resulting equations, we find that 

C + \\e\\i) 

= -- J (divv(r/ + 9) + v ■ (V77 + V9))dx -J [(U + v) ■ Vr? • rj + r/ 2 divv + ed'w(P(U + v))rj\ dx 
+ / (eF - v • Vv - 9Vr] - if79) ■ vdx + / (eG - v • V0 - (77 + 9)divv)9dx 



(2.26) 



< g IHIalMlo + 5(||v||? + ||0|| 2 ) + Q^dlPH 2 .! + ||G|| 2 i) + e||P|| 2 (||£/|| 2 + ||v|| 2 )||r;||o 
+ l|v||f + ||r ? || 2 ||^ 2 + ||v|| 2 (||^|| 2 + ||r ? || 2 )], 
where we have used integration by parts, Sobolev's inequality and the fact that 

-i J ((77 + 6>)divv + v ■ (Vr] + V9))dx = ~ J [(77 + (9)divv - (77 + 6>)divv] = 0. 



Finally, if we take 5 in (|2.26p suitably small and apply Poincare's inequality, we obtain the 
estimate (j!T25|) . □ 

2.2.2 Stokes problem 

We rewrite the momentum equations (l2.10P n as an inhomogeneous Stokes problem to derive 
the desired bounds for llvlU and || — - — -lb: 



Vt7 + V6> 

-/iAv H ? = eF - v • Vv - <9Vr/ - 77 V0 - U • Vv + CVdivv, 

e 

div v = div v, 
v = 0, on 0. 



(2.27) 



By the usual estimates for the steady Stokes problem (cf. Galdi's book j9j Chapter IV]), 
Sobolev's embedding H 2 L°° and the inequality 



we have 



v 2 + 



Vr/ + V0 



|divv|| 2 < <y||v||l + C 5 ||v|| 2 , 



<C(||eP|| + ||v • Vv|| + ||0V77||o + \\rjV6\\o 
+ \\U • Vv|| + ||divv||i). 



(2.28) 



(2.29) 
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and 



I v ll3 + 



<C7(||e#||i + ||v • Vv||i + ||eV»7||i + hV0||i 
+ ||Z7-Vv|| 1 + ||divv|| 2 ), 



which together with ([2"^5) , (j2~28j) and (f2T29|) yields 
Vr/ + V6» 



l v P + 



+ l|v||f Wo + 



i < C e (l + \\U§){e(\\F\\i + \\G\\-i) + ||v||| + \\9\\s\\vh 

1/2 



ell-Plbdl^lb + l|v|| 2 ;||'/lln 
where Cq is a positive constant. 



+ V h( l + II7/H! 



)} + ||V 2 divv|| , 



(2.30) 



2.2.3 Estimate of ||V 2 divv|| 

As in |23 } I25 |, [T2] . in order to control the term ||V 2 divv||o we divide it into the interior part 
and the part near the boundary. We remark that here we have to carefully deal with the terms 
which involve with the large parameter 1/e in (|2. 10j) . 

I. Interior estimate 

First, we derive the interior estimate of V 2 divv by using the estimate (|2.29p . Let xo be a 
Co°-function, then we have 

Lemma 2.5. There is a positive constant Cj independent of e, such that 
HlxoV 2 v|| 2 + CllxoVdivvH 2 + K\\ X oV 2 e\\l 



<c 7 



+ 



l,%l|4 



(\\U\\ 3 + llvyilr/H 2 + 6 2 (||F|| 2 + ||G|| 2 ) + e||P|| 2 (||tf|| 2 + || v 1 1 2 ) 1 1 77 1 | a + ||v|| 5 (2 . 31) 

Vr/ + V0| 2 



lill'lli 



+ ni3(||v|| 2 + ||0|| 2 ) + ||v|| 2 (||0|| 2 + ||r / || 2 ) + ||v||^ 



+ d 



Proof. We differentiate (12.101) with respect to x to get that 



v j df jr i - di(U j + v j )d j r] - tid^n - ^(r/divv) - edjdiv(P(£7 + v)), 



Wd%v k + diWd jV k - fidf ijV k - C^divv + ^ + d ^ 6 = ed .pk _ q.^q^ _ d? k v ), 



5,divv 



ediG - di(Vdjd) - d l ((rj + 0)divv). 



(2.32) 

Multiplying (f2T32l) 1 . (f2732|) 2 and (pT32"j) 3 by Xo^iV, xldiV k and xl^iO in L 2 respectively, and 
summing up the resulting equations, we find that 



nWxo^Alo + Cllxo^divvH 2 + KWxo^jOWl 

-~l Xo9 l divv(a,r ? + diO) + xldiV k d k (d iV + d^dx 



{2nxodjXod? j v k d i v k + 2(xod k XodidivvdiV k + 2KXod j Xod? j 6d i d)dx 



Xo 



di{U^ + v^)dji] + (U 3 + ff > )dj i rj + djr/divv + r/djdivv + edidiv(P(U + v)) dirjdx 



+ / xl(ediF k - di(v ■ Vv k ) - d 2 lk {hmv k dx 



+ / xbizdiG - diiVdje) - d l {{r, + e)dwv))diedx. 
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If we apply partial integrations to the above identity, employ Sobolev's and Young's inequalities 
and the fact that 



Xo<9;divv(<9j77 + + Xo^iV 'd k (dir] + diO)dx 



1 



2xotoc^ fc (c> 4 77 + 9i9)dx + - / x 2 odiV k d k (d iV + 80) - x&v* 8^ + 8 l 6)dx 



+ c 5 \\Al 



Vr/ + ve 



= ^ J 2xod k XodiV k (dir] + 8i9)dx < 5 
we infer by summing up i,j, k that 

^I|X0V 2 V||2 + C || X()Vdivv ||2 + K || Xo V 2 0||2 



<OTIs + l|v||3)IHI? + « 



+C S 



+ 



11 + 







i + l|v||f + 



^(11^115 + lW + e||P|| 2 (||C/|| 2 + ||v|| 2j 

+ Il^ll3(||v||? + ||0||?) + ||v||2(||0|| 2 + || ?? || 2 ) + ||v||?], 

which, by using Poincare's inequality and choosing 5 appropriately small, implies the lemma. □ 
Lemma 2.6. There is a positive constant Cs independent of e, such that 
HlxoV 3 v|| 2 + C||xoV 2 divv|| 2 + K\\xoV 3 e\\ 2 Q 

3 + ||v|| 3 )|k||! + t 2 (\\F\\l + IIGH 2 ) + e\\PH\\U\\ 3 + ||v||3)|k|| 2 + ||v||S 

Vt/ + V0 2 



<C 8l 

+ MM + WHhg + + ||v||!j(||0||jj + llr/lH) + ||v ll2 
Proof. We differentiate (|2.10p twice with respect to x to get that 



(2.32) 



+ 6 



u j d? ljV + 



<9,divv 



- dftiW + &)d jV - 8(W + v>)dj lV 
-di{W + v j )d? jV - ^/(r/divv) - e8ldw(P(U + v)), 



W8f jlV k + 8iWdf jV k + d\W8jV k + diWdfjV* - ^df ljjV k - C^ fc divv + 
^ -df^d^-dUn), 



dfikV + df lk e 



ed?,F k 



(2.33) 



o2 j • 

wdf^e + die + dlwdje + diWdf.e - k8% j + 



edi l G-df l (Vd 3 e)-dt l (( v + e)dw^. 



Multiplying (plgg]^ . (|2~33"1) 2 and (pOBl o again by Xo<9 2 ^, Xo^* and Xo<^6> respectively, and 
summing up the resulting equations, we deduce that 

MllXo^ll 2 , + Cllxo^divvH 2 + 4xody\\ 2 
= ~\f xldudi^dan + dle) + xld 2 a v k d k {d 2 u n + dl9)dx 

2 X od 3 x^dl jV k dlv k + Kdlfldle) + 2(xod kXo d u divvd u v k dx 

Wdf ljV + vWl jV + d\ (W + vi)d jV + 9,(17' + vi)d] lV 

+8 l (W + vi)df jV + didivv) + ed 2 div(P(C7 + v))] Xpfadx 
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+dlWd jV k + diWdf jV k ) 



+ I ed\F k - a&tfdjv*) - dl k {h) ~ (U^df jlV k + d t Wdf jV k 

xld 2 u v k dx 

ed\G - dl^djO) - 4(( V + 0~)divv) - (Wdf^e + d x Wd%e 

+dlU j djd + diWdfjO)] xlduOdx. 
We integrate by parts the above identity, utilize Sobolev's inequality and the fact that 

X^divv^r? + die) + xldlv k d k (dlv + dl6)dx 

\ I xldf l v k d k (dl7 1 + dl6)-xldlv k d k (dlv + df l 9)dx + - e J 2 X od kX odlv k (df lV + d 2 e)dx 

Vt] + V6 2 



- / 2x d k xodlv k {&( l 7 1 + &i l e)dx<5 



+ C S \\v\\i, 



and sum up k to infer that 



^HxoV 3 v|lo + CI|xoV 2 divv||2 + K || XoV 3 0||2 



< (11^113 + [|V|| 3 )IN1+ «(||V||I + 



l§+ 



Vr/ + V0 



+C s [e z (\\F\\i + \\G\\i) + e||P|| 3 (||tr||3 + ||v|| 3 )||r/||2 + Mt + WvMPM 



+\WH\w\\i + Pm + \\nimr2 + \\riM) + hr2\, 

which, by employing Poincare's inequality and choosing 5 suitably small, gives the lemma. □ 
II. Boundary estimate 

Next, we shall use the method of local coordinates to bound V 2 divv in the vicinity of the 
boundary (also see [23j [23 E2])- For completeness, we briefly describe the local coordinates 
as follows. First, one construct the local coordinates by the isothermal coordinates X((p, <j>) to 
derive an estimate near the boundary (see also |23[ [25] ) , where 



A^ • Xip > 0, • A^, > 0, Xip ■ \q 



0. 



The boundary d£l can be covered by a finite number of bounded open sets W k C M 3 , k 
1, 2, • • • , L, such that for any x £ W k D £1, 



x = A k (<p, <f>, r) = A fc (y?, 4>) + m(X k (ip, </>)), 



(2.34) 



where X k (f, 4>) is the isothermal coordinates and n is the unit outer normal to dQ. 

Without confusion, we will omit the superscript k in each W k in the following. We construct 
the orthonormal system corresponding to the local coordinates by 



ei := 



A 

I A, 



&2 



I A„ 



e 3 := ei x e 2 = n(A). 



(2.35) 



By a straightforward calculation, we see that for sufficiently small r and J £ C 2 , 

dx 



J := det JacA = det- 



A^xAj,' e 3 



J d(</?,c/>,r) ' v 

|A V ||A^| + rdA^ln^ • e 2 + jA^n^ • e%) + r 2 [(n^ • ei)(n^ • e 2 ) - (n^ • e 2 )(n (j) ■ e%)] > 0. 
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And, we can easily derive the following relations as (JacA ) o A = (JacA) (also see |23j): 

[V(A- 1 )>A= J-^A^xea), (2.36) 
[V(A- 1 ) 2 ]oA = J- 1 ( e3 xA,) 1 (2.37) 
[VtA- 1 ) 3 ] o A = J-\A V x A ) = e 3 , (2.38) 

where the symbol o stands for the composite of operators. Set y := (ip,(j),r), and denote by 
Di the partial derivative with respect to yi in local coordinates. We set the unknowns in local 
coordinates 

fj(t, y) := V (t, A(y)), v(t, y) := v(t, A(y)), y) := 0(t, A(y)), 

and the knowns 

U(t, y) := U(t, A(y)), $(t, y) := $(t, A(y)), l(t, y) := 0(f, A(y)). 

Then, we rewrite the system (|2.10|) in [0,T] x 0, where := A~ l (W n O), as follows. 

^ » afc -OfcU- 7 - • - • ..... 

U j a kj D k fj H ^- = -v 3 a kj D k f) - fja kj D k v J - ea kj D k (P(U J + £ J )), 

U 3 a kj D k v l - fiakjDkiaijDiv 1 ) - CakiD^aijDiv 3 ) + afa - Pfc ^ + 6> ) 
= e.P - v j a kj D k v l - 6a ki D k i] - f)a ki D k 6, 

U j a kj D k 9 - Ka kj D k (aijDi6) H ^— = eG - v J a kj D k 8 - fja kj D k v J - Oa k jD k v> , 

with boundary conditions 

v(t,j/) = 0, §(t,y) = on 90, (2.40) 

where a^- is the (i, j)-th entry of the matrix Jac(A _1 ) = Clearly, a^- is a C 2 -function, and 
it follows from (|236]) - (|2T38D that 

3 3 3 

a3j«3j = |n| 2 = 1, oij03i = a2i«3i = 0. (2.41) 
j=i j=i j=i 

Moreover, this localized system has the following properties (see also |23j): 

Proposition 2.1. Di(Jciij) = 0, for j = 1, 2, 3; <^D T v = 0, <;D T D^ir = on 90 in f/ie tangential 
directions r,£ = 1,2, w/iere ? 6 C^°(A _1 (Ty)). Similarly, qD T 6 = 0, qD T D{.6 = on 90. 

Recalling = X^=i a j«^> we W1U frequently make use of the following relations without 
pointing out explicitly in subsequent calculations: 

HA/v|kp(n) < C\\^xv\\lp(u), \\Dyv\\ LP (n) < C\\V x v\\ w i,p^, 1 < p < oo. (2.42) 

The above inequalities apply to rj, 6 and f7, v, 9, too. 

By virtue of the interpolation || • \\ 2 H2 < <5|| • ||^- 3 + Cg\\ ■ \\ 2 H i, the boundary estimate of 
||V 2 divu|| L 2( L 2) can be reduced to the boundedness of 



f I Jx^DKajiDjU^dyds, 
JO JQ 



where x 1S a Co°(A 1 (VF))-function. So, we can split the estimate of derivatives on the bound- 
ary into two parts: the estimate of derivatives in the tangential directions and in the normal 
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direction. 

Part 1. Estimate of derivatives in the tangential directions 

First, we apply D 2 ^ to (|2,39p with r, £ being the tangential directions to d(l to get 

' - 1 
UJa kj D lr^ + -D 2 TS [a kj D k V] 

= D^Wa kj )D 2 T fj + D^(Wa kj )D k f) + D T (W a kj )D 2 k( f, 
-D 2 ^a kj D k fj + T]a kj D k vi + ea kj D k (P(W + &))], 

Wa kj Dl T ^ - D^a^D^D^) + (a ki D k ( aij £>,#')] + + <?)] 

+D^[e-F l - v j a kj D k v l - 9a ki D k i) - rja ki D k 9], 

Wa kj D 3 k ^6 - KD^ayDkicujDfl)] + ±D^[a kj D k v>] 
= D^Wa kj )D 2 J + D 2 ^(Wa k3 )D k § + D T (Wa kj )D 2 ^ 

+D 2 Tl \eG - v j a kj D k 9 - i)a kj D k v j - 9a kj D k v j \. 

We multiply ^35D i, d233D 2 and (g3SD 3 by J X 2 D T ^fj, J X 2 D Tf y l and J X 2 D ri :9 respectively, and 
integrate the resulting identities to deduce that 

D^akjDkiaijDiv 1 ) + (akiDkiatjD^)] • J X 2 D 2 l .v i dy 

- [ KDlAa^DkiatjDiO)] ■ J X 2 D 2 Jdy 
Jn 

+ j_ Wa kj (Dl^ ■ J X 2 D%f, + Dl^ ■ J X 2 D%v i + D\J ■ J X 2 D 2 J)dy 

+- / [DL[a kj D k v>] ■ Jx 2 D 2 T( .f, + [ aki D k (rj + (9)] • Jx 2 L>^ + D 2 Aa kj D k v^} • J X 2 D 2 J}dy 

e ./fi 

= ^{^(^a fcj )Z)L^ + ^(^a^)^^ + ^r(^a fei )^r} (2.44) 

-D 2 T ^a kj D k fi + rja kj D k v* + ea kj D k (P(W + #))]} • J X 2 B\f]dy 

+ I {D^a^Dl^ + I%Jfra kj )D k v i + D T (Uia kj )I&tr i 
Jn 

+D 2 .^eF i - i j a kj D k tf - 9a ki D k fj - rja k iD k §}} • J X 2 D^v i dy 

+ I {D^Wa kj )D 2 k J + D 2 T JWa kj )D k 9 + D T (Wa kj )D 2 k J 
Jn 

+D 2 ^[eG - &a kj D k 9 - fja kj D k & - §a kj D k &}} ■ J X 2 D 2 r( 9dy. 

Now, we denote LHS of (|2.44p :=Z4 + L2 + £3 + £4 and have to deal with each term due to 
integration by part and the boundary conditions. 

+D^(fia kj )D k [D T (ai j )Div l + ayD^u*] + fJ,a kj D k ^ [D T {aij)Div i + ayZ&u*] 
+D 2 ^Cau)D k (a lj D l n>) + D T ^a ki )Dl i (a lj D l v i ) 
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2 ,V'l 



where 



and 



+I> c (Co w )I> fc [D T (ay)Av J + aijDfr 
+Ca ki Dlz[D T (a lj )D l v> + ayAr^']} ■ Jx 2 D 2 T ^dy 

- f \ f M kj Dl i (ai j Df T v i ) + Ca ki D 2 H ( aij Dl^)} ■ J x 2 D 2 ^v i dy 
Jn 

fiJx 2 a kj Dl T ^v i a lj Df T i:V i + CJx 2 akiDl rl .v i a lj D^v j dy 

+ / [D k (fiJx 2 a kj )D i (ai j D lT v i ) ■ D^v* + fiJx 2 a kj D 3 kr ,v i D^a lj )Df T v i 
Jn 

+D k (QJ X 2 a ki )D i (a lj D 2 lT v i ) ■ D 2 .^ + fiJ^akiD^^D^ai^Df^dy, 

L' 2 = Kj X 2 a kj Dl T Ja lj D? T Jdy 
Jn 

+ / \D k {KJx 2 a k] )D^a l3 Dld)D 2 T ^e + K^a^Dl^D^a^Dfj] dy 
Jn L J 

D T^akj)D k (aijDi§) + D T {Ka kj )D 2 H {a l:j Di9) 
+ D^na kj )D k [£> T (ay)A0 + ayl&fl] + Ka kj D 2 K D T ( aij ) D t §} ■ J X 2 D 2 ^dy. 
On the other hand, recalling that a k jD k \ji = 0, we have 

4 = - \ f Jx 2 a kj D k W(\D 2 ^\ 2 + iD^fif + \D Ti 6\ 2 )dy 

= -\i D k (J X 2 a kj )W(\D 2 ^\ 2 + \D%vf + |D?^| 2 )dy. 

As for L4, in view of the following identity 

- f D 2 ^[a ki D k {f, + §)] ■ J X 2 D 2 T( v i dy = -- f J X 2 D 2 ( (a kt D k v i )D 2 ( (rj + 9)dy 
e Jn e Jn 

D 2 T( (a ki )D k (f, + §) + D T (a ki )D 2 z(f) + 9) + D^a ki )D 2 T (fj + 0)] • J X 2 D 2 ^dy 

6 7n 



+- 

e Jn 

1 



-D k (J X z )a kl D^v l - D k (a ki )J X z D z Ti :V l ]dy, 
we deduce that 



4 =- / \D 2 Ti (a ki )D k (f} + §) + D T (a ki )D H (fj + 9) + D^a ki )D 2 T (fj + (9)1 • J X 2 D%v l dy 
e Jn L J 

+ 1 jT £)| T W + ^[^(afciJDfcfi* + Dria^D 2 ^ + D^D 2 ^ 



- D k (J X 2 )a ki D^v l - D k {a kl )Jx z D^v l ]dy. 
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Substituting the above estimates into (|2.44j) . using Sobolev's and Young's inequalities and 
taking into account the property (|2.42[) . we deduce that 



nJX l a kj Di r zV i a lj D$ ri v t dy + / CJ^a hi D% Te v g a lj Df Te v 3 dy + / KJx^ahjD^OaijDf^edy 



•>3 -i. 



<C 9 



■)3 ~i 



3 



:i!K/l|2 + Il v ll3 



2 , " -»| -I- e[|^[| 3 Cll^ll3 -H l|v[| 3 



2 + 

2n\~\\2 , imi|2\ , n~i|2 



l(l|v||i + 



+ e z (\\F\\t + \\G\\i) + wnt + wvimm + mi) + wm 



+ v ? + 



+ 6 \W\\i + 



11+ 



(2.45) 



where Cg is a constant. 

Part 2. Estimate of derivatives in the normal direction 
We multiply (|2.39|) 2 by a 3 i to obtain that 

- (p+ QDz{a l3 Div j ) + ^D 3 {f) + 9) 

= - a 3i U j a kj D k v l + ea 3i F - a 3i v j a kj D k v - 9D 3 t) - t)D 3 6 
+ V a 3i a kj D k (ai j Div t ) - D^ayDfl) 
where the last term in RHS of (|2,46p can be written as follows, 
/i a 3 ia kj D k (aijDiv % ) - D^ayD^) 

= n D 3 (a 3j )D 3 v j + D 3 {a Tj )D T v j + a Tj Dl T v j - a 3j D 3 {a 3 j)a 3i D 3 v l 



(2.46) 



-a T ja 3i D T aijDiv l - a T ja^a 3i D T ^v l - a 3 ja 3i D 3 (a T j)D T v 

which does not include the term D 33 v. 

Step 1. To continue our estimate, we show the following lemma. 

Lemma 2.7. There are a constant C\q and a small 5 > 0, such that 



r, £ = 1,2, (2.47) 



2 

<C W { 
+ 



= Jx 2 \D z T3 (a lj D l V)\ z d y + - I J X 2 \Di 3 (a kj D k e)\'dy 



v||? + e 2 \\F\\'i + ||v||4 + \\0\\l\\ V f 2 + / J X 2 \D* £v v| a dy 



+ V 3 



+ + 



.iJ||<7||l + €||P|| 3 (||Cr||3 + ||v||3 



(2.48) 



+ e 2 ||G||? + ||v| 



+ v 



+ 5 v i + 



ID- 



Proof. We differentiate (I2.46P with respect to y T (r = 1,2), then multiply —Jx 2 F> r3 (aijDiv^) 
in L 2 (fl) to get 



J X z \DUaijDiv j )\ 2 dy - - I Jx'DWn + e)D^(a lj D l v^)dy 



e Jn 



< 



C(\\a 3i Wa kj D k v l \\i + e z \\F l \\i + \\a 3i Va kj D k v 
+ \\§D 3 fi + t?L> 3 I||?) +C f J X 2 \D% y v\ 2 dy 



2 
li 



(2.49) 



<C(MlMl5 + *W + W + W\i\\rj\\i) + C I J X 2 \D^ y v\'dy 
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(2.50) 



In the mean while, we apply D t3 to (|2.39P o and ([2.39^ . take the product of the resulting 
equations with J X 2 D 2 3 9 and J X 2 D 2 3 i) in L 2 (Cl), and sum then two identities to get 

- f kD t3 \a kj D k {flijDiB)\ ■ J X 2 D 2 3 9dy + - [ D 2 r3 (a kj D k v j ) ■ J X 2 (D 2 r3 9 + D 2 r3 fj)dy 
Jn L J e Jn 

+ / [D 2 3 (Wa kj D k ri) ■ J X 2 D 2 T3 f, + D 2 3 (Wa kj D k 9) ■ J X 2 D 2 3 9]dy 
Jn 

= - J D 2 T3 \i j a kj D k f, + fia kj D k d j + ea kj D k P(U j + v j ) } • J X 2 D 2 

+ j D 2 T3 (e& - d j a kj D k 9 - f]a kj D k d j - 9a kj D k ^ ■ J X 2 D 2 3 9dy. 

We denote LHS of (|2.50p := L'[ + L'^ + L 3 . To control L' k , we integrate by part to deduce that 

Li" =k [ J^DlsiakjDk^Dl^D&dy 
Jn 

KJ^D^e^Dl^a^Dkia^O) + D T (a k j)D 3k (aijDi9) + D 3 (a kj )Dl T (a lj Di9)Jdy 

+ k [ DkiJx^D^D^aijD^dy 
Jn 

- k f^J X 2 [D 2 T3 {a k3 )D k 6 + D T (a kj )D 2 3k 9 + D 3 (a kj )D 2 T 9^j ■ D 2 3 ( aij Di§)dy 



and 



W = f_ J X 2 D 2 3 v ■ [D 2 r3 {Wa k] )D k f, + D T (Wa kj )D 2 3 fj + D 3 (W >a kj )D 2 T i))dy 

~ \ ! D k (J X 2 a kj )W\D 2 r3 fi\ 2 + J X 2 a kj D k W\D 2 3 \ 2 dy 
1 Jn 

+ J (D 2 T3 (U^a kl )D k 9 + D T (Wa kj )D 2 k3 9 + D 3 (W a kj )D 2 T 9) ■ J X 2 D 2 r3 9dy 

- \ ! D k (J X 2 a kj )W\D 2 T3 9\ 2 + J X 2 a kj D k W\D 2 3 §\ 2 dy. 
1 Jn 

Inserting the estimates for L'{ and L 3 into (|2,50p . we obtain 

£ / J X 2 \D 2 T3 (a kj D k 9)\ 2 dy+ - f J X 2 : D 2 3 tf + 9)D 2 T3 { ai] D^)dy 
1 Jn e Jn 



< Ml + 5[ WDKa^fo + \\D 2 T {aijDi9)\\ 2 + \\9\\ 2 



+C 



3\\m\2 + \\uh\\e\\^ + \M 3 \\vU + e\\Phm\3 + \\yh 



+e 2 \\G\\i + 11^11111^111 + ll^llIlMiJ- (2-51) 

Thanks to Sobolev's and Young's inequalities, we take the sum of (|2.49p and f|2.51j) to deduce 
the estimate (EOHD . □ 

Step 2. Now, it suffices to bound \\D 33 {aijDiV^)\\Q in order to close the estimate for divv. 
We apply D 3 to fl!P6D to find that 



19 



= D 3 (a 3i )(-U j a^D^ + eF i - v^ojy - a 3i D 3 (U j a kj D k ^ ~ + Va^D^) 
-D 3 0D 3 fj + f)D 3 6) + 0(l)(Dl 3T & + D 2 ^' + D&). (2.52) 
Now, multiplying the above equality (12 .521) by —Jx 2 D 33 (aijDivi) in L 2 (f2), one infers that 



J^l^isK-A^')! 2 ^ - - / Jx 2 DUa lj D l iP) ■ D 2 33 {f, + 9)dy 



<C(ll^ll5lHl5 + ^||F||f+||t;||5 + 



llll^lll) + IMI1 



HI 



(2.53) 



J X 2 \D 3 33r ^dy. 



Correspondingly, applying D 33 to (|2.39[) i and ()2.39j) ^ and multiplying the resulting equations 
by Jx 2 D 33 fj and Jx 2 D 33 fj respectively, we get 



_ J X 2 \D 2 3 (a kj D k 6)\ 2 dy + - J _ J X 2 D 2 33 ( aij D^) ■ D 2 3 (fj + 6)dy 



(2.54) 



<C\\9\\i(l + \\U\\ 3 ) + 6\\9\\i + \\U\\ 3 \\r,\\i + \\v\\ 3 \\ V \\ 2 2 + e||P|| 3 (||t/|| 3 + \\v\\ 3 ) 
+ e 2 ||G||f+||S|||||0|H+||S||I|H|§. 

Combining (]2.53p with (|2.54p . we see that there are a constant C\\ and a small S, such that 



■_ Jx 2 \D 2 33 (a lj D l V)\*dy + ^ I J X 2 \D 2 33 {a kj D k 9)\ 2 dy 



,2 1 t->2 



< 



2 -4- 

3 + 



|v|||) + Cn{(||^||i||v|| 2 + e 2 ||F|| 2 + ||v||^ + ---^ 



+ 



v||?+ / J X 2 | J D 3 3 3.v|^y + ||0||<f(l + 



3) + 



+e||P||3(||CA||3 + ||v| 



2 + e 2 ||G|| 2 + ||v|| 2 ||^|| 2 + ||v|| 2 ||r ? || 2 



3lklll + l|v|| 3 | 

.Ml 2 i . . i i 12~| I 



(2.55) 



Step 3. To control the term -D| 3r v on RHS of (12351) . we introduce an auxiliary Stokes 
problem in the original coordinates in the region near the boundary: 

-fj,A x [( X D T v) o A- 1 ] +^V x [(xD T (v + e))oA- 1 } =G t in Wntt, 

dw x [( X D T v) o A- 1 } = G 2 in WnQ, 
(xA-v) o A -1 = in Wnfi, 



where 



G\ =xD r Ca ki D k (aijDiv j ) + eF l -~v j a kj D k v l -~9a ki D k f) - fja ki D k 9 



U^a kj D k v l + o(l) Div % + D z kl v l + -D k (f) + 



and 



G l 2 = o(l)(D T V + D k V + D 2 Tk V) 
which can be bounded as follows. 

V(v + 9) 2 



HGi||| 2 <C( 



L 2 



+ I J X 2 \D 2 Tk (a lj D l vi)\ 2 dy) + 5\\v\\ 2 



(2.56) 



C 5 ||v|| 2 + C( e 2 ||F|| 2 + ||v||l + ||0|| 2 ||r?|| 2 + ||L^||v|| 2 ) 



and 



\\G 2 \\ 2 H i < *||v||l + Q||v|| 2 +C f Jx 2 \D Tk ( aij D^)\ 2 dy. (2.57) 

Jn 
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Due to the regularity theory of the Stokes problem (see [U] ) , one has 
/ \A x ( x D T v) o K-\ x )\ 2 dx < CCIIdHi^nn) + ||G 2 || 
where the left-hand side of (|2.58j) is equal to 

p a 3 3 3 

/ \/\ x ( X D T v)ok-\x)\ 2 dx= / j| V Va fci £ fe (yVA(x£>rV)) 

3 

= / E a fci %^ T v| 2 d2/ + (l) [(\D T v\ 2 + \D 2 yr v 

n j,k,l=i n 

And we use (j'2.4ip to get 



(2.58) 



dy 



l )dy. 



3 3 

^33rV= Yl a kj a lj D llrV- ^ E a kj a lj D llr^^ 
j,k,l=l l<k,l<2j=l 

from which, (12.56P and (12.57P it follows that the inequality (12.58P gives 

(Cu + 1) / Jx^fsrVl 2 ^ 

Jn 

Cdl^ill^^nn) + \\G2\\ 2 HHwnn) ) + C f J X 2 \D% c \ 2 dy + C 5 ||Vv|| 2 2 + 5||v 

V(t? + 0) 2 



< 



|v||? + 



<5\W\\ 2 3 + C 12 

+ (e 2 ||^lll + l|v||^ + ||^||2ll??ll2 + 

Now, letting 



e 

|2 1 1 II 2 



+ / Jx 2 (|£>f r ,v| 2 + \Dl T { ai iDiV)?)dy 
Ja 



4 1| || 2- 
3 ll v lll, 



(2.59) 



and 



*x 



J X 2 (\D% y v\ 2 + \D 2 3 ( aij D^)\ 2 + |^| 3 (ayA^)| 2 + l^ 3 3 3rv| 2 )^ 
J X 2 (a kj \Dl T ^\ 2 + |D 2 3 (a fcj D fc 0)| 2 + (^(a^D^)! 2 )^, 



we can apply Cauchy-Schwarz's and Young's inequalities as well as the estimate (|2.29|) to deduce 
from (H35|), (H3HD, (12351) and (1239]) that 



alHli + (l|v|| 3 + HvllDhUi + e||P|| 3 (||C/||3 + ||v|| 3 j||v/i > 



+ C(1 + 



!)(l|v|| 2 + 



) + S ||v||| + 



11 + 



V(r/ + fl) 2 

e l 



(2.60) 



+ C 5 



e 2 (ii^ii 2 + iigii?) + \\m + mi + iivii 2 )hn 2 + wnin! 



which, together with (]2.30p - (j2.32|) . results in 



|v|ii + 



11 + 



Vri + V6 2 
1 



l + l|v||| + ||v|| 3 ) 



<C 12 {||C/|| 3 |MI + (1 + \\U\\I) [e||P|| 3 (||C/|| 3 + ||v|| 8 ) + ( 

+ (1 + ll^ll!)(l|v|| 2 + + (1 + \\U\\i) [e 2 (\\F\\ 2 + ||G|| 2 ) + + ||v|| 2 ||0|| 2 ] } 



(2.61) 
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2.2.4 Boundedness of rj 

In the next lemma, we derive upper bounds of \\r]\\i and ||??|| 2 . 

Lemma 2.8. There are a small 5 > 0, and two positive constants C13 and C14 independent of 
e, such that 



<Ci 3 



e 2 ||F|| + e (||v|| 2 ||v||i + Whhh + ni2||v|| 2 + ||diw||i + (||F||-i + ||G||_i)) 



+ e || j P|| 2 (||^|| 2 + ||v|| 2 ) + iiviniMio + ||v||5 + \\r,U\Qh + ||v||i(||0||i + N|i) (2.62) 



and 



< C 



+ 



u (l + e)(l + \\Ug){e(\\F\\ 1 + \\Gl 



■i) + l|v||i + 



2+ V 



e||P|| 2 (||C/|| 2 + ||v| 



2; '/ 



+ v 1 



i + Wl)}+e||v|| 3 + <5||0|| 3 . (2.63) 



Proof. Prom (|2.29p we get 

e||v|| 2 + ||Vr/ + V9\\o < C[e 2 \\F\\ + e (||v|| 2 ||v||i + H^alhlb + ||^|| 2 ||v|| 2 + ||diw||i)], 

which together with Lemma 12.41 gives 

[|V77||i<||V»7 + Ve||o + [|Ve[|o 
< C[e 2 ||F|| + edlvbUvIK + \\§\\ 2 \\ V \\ 2 + ||C/|| 2 ||v|| 2 + \\dwv\U)] 



+C 



vblMI 2 , + e 2 (||F|| 2 1 + \\G\\ 2 _ t ) + 6||P|| a (||l7|| 2 + ||v|| 2 

1/2 



+ l|v||f + llr/^ll^ + llvllfdl^^ + llr?!! 2 
If we apply Poincare's and Young's inequalities to the above inequality, and use the fact that 
{A l +A 2 + ---+A n ) 1 / 2 < A{ /2 + A] /2 + ■ ■ ■ + A}/ 2 for^>0 (t = 1 



(2.64) 



we obtain the estimate (|2.62j) immediately. 

On the other hand, from the estimate ()2.30p we conclude that 



\\Vvh < \\Vv + ve\\i + \\wh 



<C6(1 + ||C/|| 2 ){6(||F|| 1 + ||G||_ 



•i) + l|v||§ + 



2 + V 



+ 



e||P|| 2 (||[/|| 2 + ||v|| 



+ l|v|| 1 (||e||i + |H|i)}+e||v|| 3 + *||6>||3 + C',||e|| 1 , 



which, together with Poincare's inequality, (|2.62p and (|2.64p . implies (|2.63p . 



□ 



3 Existence of the nonlinear problem 

In this section, we give the proof of the existence for the nonlinear problem (|1.5p by using 
the Tikhonov Theorem which can be found in [20]. For completeness, we state the theorem in 
the following. 

Theorem 3.1. (Tikhonov Theorem, [20, P72, 1.2.11.6]) Let M be a nonempty bounded closed 
convex subset of a separable reflexive Banach space X and let F : M — )• M be a weakly continuous 
mapping (i.e., if x n £ M , x n — 1 x weakly in X, then F(x n ) — F{x) weakly in X as well). Then 
F has at least one fixed point in M . 
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Define a Banach space X by 



X = H 1 x Hq x Hi, 



which can be easily verified to be separable and reflexive. 
A convex subset K\ (E) of X is defined by 

Ki{E) = | (v, 9) E (P 3 n Hq)) x (H 3 n Hq) I ||v|| 3 + ||0||3<P}, 

where E < 1 is a small positive constant. By the lower semi-continuity of norms, we easily see 
that the subset K\{E) is also closed in X. 
We define a space K by 

K = K xK 1 (E), 



where Kq is defined by f|2 . 9j> . Note that K is a nonempty bounded closed convex subset of X. 
Now, we define a nonlinear operator N from K to X by 

N{U,v,9) := (U,v,0), 

where U and (v, 9) are the solutions of (|2.3p and (|2.10p for given (U,v,9), respectively. 

Next, we want to find a fixed point (U,v,9) of iV in K, such that (U,v,9) = N(U,v,9), 
which, together with the existence of weak solution in Lemmas 12.11 and 12.41 gives that (U, P) 
and (n,v,9) are solutions of the boundary value problems (|2.ip and (|2.2p . respectively. So 
(U + v, eP + 7], 9) will be a solution to (II. 5h . For this purpose, we have to show that N maps K 
into itself and N : K -4 K is a weakly continuous mapping. 

Lemma 3.1. There is a small constant eo > 0, depending only on f2,^,A, f and g, such that 
for any e G (0, eo), K is a nonempty bounded closed convex subset of X and N(K) C K. 

Proof. By virtue of the definition, it is obvious that K C X is a nonempty, bounded, closed 
convex set. Now, we will show that the operator ./V maps K into itself, i.e., N(K) C K. To this 
end, let (U,v,9) C K and (U,v,9) = N(U,v,9). By Lemmas O and we see that U <E K 
for all [> G K . Thus, it suffices to check that (v,0) G ^i(^) for (v,6>) G i^i(-E). By (|27FD and 
(HS|), we have 

+ ||VP||i < Ml and ||C/|| 4 + ||VP|| 2 < M 2 , (3.1) 

where Mi = C 3 ||h||i(||h||i + l) 8 and M 2 = C4||h|| 2 (||h|| 2 + l) 12 . 

On the other hand, recalling the definition of F and G, we get from (I3.ip that 



II^Hi =|| (eP + r/)f - (eP + ??)([/ + v) • V(U + v) - 6>VP - PW||i 

<C[|M| 2 (||f ||i + ||LT||1 + ||v|| 2 ) + e||P|| 2 (||f ||i + \\U\\l + ||v|| 2 + ||flf|| 2 )] (3.2) 
<C[|r7|| 2 (||f ||i + (Mi + l) 2 ) + eCMi(||f||i + (Mi + l) 2 ), 

||G||i =||* - (eP + rj)(U + v) • V9 + (eP + r/)0divv + Pdivv||i 

<C[e(\\U\\ 2 2 + l|v|||) + (h|| 2 + e||P|| 2 )(||t/||2 + ||v||a)||e|| 2 + l|P||2||v|| 2 ] (3.3) 
<C||r/|| 2 (Mi + 1) + C(e(Mi + l) 2 + (M x + 1)). 
As a result of Poincare's inequality and Lemma 12.41 we have 

||v||i + ||0||i < C[(E^ + £7)||77|| 2 + e[|M| 2 (Mi + 1) + C(||7/|| 2 (||f ||i + {M x + if) 

+eCMi(||f ||i + (Mi + l) 2 ) + e(Mi + l) 2 + (Mi + 1)) 

+eC 5 (Mi + l) 2 + £ 2 } + 5\\n\\ 2 
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< C[E\ + e(||f lli + (Mi + l) 2 )] \\ V \\ 2 + e 2 CMi(||f ||i + (Mi + l) 2 ) + CE 2 
+5\\7 ] \\ 2 + eC 5 (M 1 + l) 2 , (3.4) 



where we have used the estimates (|3.ip . (|3.2p and f|3. 3|) . 

On the other hand, in view of Poincare's and Young's inequalities, (|3.4p and f)2.63j) in Lemma 
we find that 

1Mb < C(l + e)(l + Af- 1 2 ){e[||7 7 || 2 (||f|| 1 + (Mi + l) 2 ) + Mi(||f \U + (M 1 + l) 2 ) 
+h\\ 2 {M x + 1) + (e(Mi + l) 2 + (Mi + 1))] + E 2 + (E + E/5)||t7|| 2 
+e^Mf (Mi + £7)3 II^HI + £(£ + H77II2)} + e||v|| 3 + <5||0|| 3 
< Cib(1 + e)(l + M 2 ){e[||f||i + (M x + l) 2 + (M x + 1)] + (£7 + £75) + 5}||r?|| 2 
+e||v|| 3 + a||0|| 3 + Ci 5 (l + e)(l + M 2 ){e(Mi + l) 2 



Mi(||f ||i + (Mi + l) 2 ) + (e(Mi + l) 2 + (Mi + 1))] + £7 2 }, 



(3.5) 



where C15 is a positive constant. 

Combining (|2.6ip with (j3.2|) - (j3.5|) . we conclude that there is a constant C\%, such that 

IMIs + ||0||s + hh < Cie(l + Mi) 5 {e [||f ||i + (Mi + l) 2 ] + (E + £75) + <s}|H| a 

+eM^||v||3 + ^||0||3 + Ci6(l + Mi) 5 {e(Mi + l) 2 



+ £7' 



(3.6) 



+e(l + Mi)[||f||i + (Mi + l) 2 

Thus, first taking 5 small enough and then choosing eo and E suitably small, such that 

Ci 6 (l + Mi) 5 {e [||f ||i + (Mi + l) 2 ] + (£7 + £75)} < 1, e Mf < 1, 

Ci 6 (l + Mi) 5 {e (Mi + l) 2 + e (l + Mi) [||f ||i + (Mi + l) 2 ] + £ 2 } < E, 

we deduce from (|3.6p that for all e G (0, eo), 

IMI3 + ||0[| 3 + M 2 <£, 

which gives 1 1 v 1 1 3 + ||0||3 < E immediately. This completes the proof. □ 

Lemma 3.2. Let N,X,Kq and K\{E) be the same as in Lemma \3.1\ Then N : K — > K is a 
weakly continuous mapping. 

Proof. By the definition of weakly continuous mapping (see, for example, [201 P72,1.4.11.6]), it 
suffices to prove that N is continuous on K in the norm of X. 

Let (Ui,Vi,9i) = N(Ui,Vi,9i), i = 1,2. In particular, let (Ui,Pi) G (i/ 4 n H^) x H 3 and 
(Vi, Vi, di) e H 2 x (H 3 nff(J)x (£T 3 n F 1 ) be the solutions of ([23]) and (|2TTUj) for given (t/i, v i} 0*) 
respectively, i.e., 

I (tT< + vO • VUi - nAUi + VP; = h, J Pdx = 0, 
\ div Ui = 0; 



(3.7) 
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and 

Ui ■ Vr/i + ^p- = -Vj • Vr/j - ??idivvi - ediv(Pj(£/j + Vj)), 

Ui • Vv, - M A Vi - CVdivv, + V??t + V6i = eF t - y t • Vv, - e^m - vNOi, ( 3 - 8 ) 

e 

C/j • VOi - nA9i H - = eGi - v» • V0, - r/idivv; - ^divVj, 

e 

where the force and heat source G{ are given by 

Fi = (ePi + TH)f - (ePi + mWi + v t ) ■ V(Z7< + v*) - ON Pi - PW,, 
Gi = ^i- (ePi + Vi )(Ui + Vj) • V6i + (ePj + r/^divv; + Pdivv;. 

Now, if we set 

w = u 2 -u 1 , w = u 2 -u u q = p 2 -p 1 , t = m-m, 

w = v 2 -vi, w = v 2 -vi, (3 = 2 — 0i, P = 02 — 01, 
J = F 2 — F\ , I = G 2 — Gi, 

then, we can have the following systems: 



{Ui + vi) ■ VW - fiAW + VQ = -(W + w) • W 2 , J Qdx = 0, 
divW = 0, 



and 



(3.9) 



( divw 

Ui ■ V£ + = -div(vi^ + wr] 2 ) - ediv(Pi(VF + w) + Q{U 2 + v 2 )) - W • Vt? 2 , 

V£ + V/3 

Z7i • Vw - //Aw - CVdivw H - 

= eJ - w • Vv 2 - vi • Vw - W • Vv 2 - V(0i£ + /3r/ 2 ), (3-10) 

Z7i • V/3 - kA/3 + 

_ e _ 
= e/ — w • W 2 — Vi • V/3 — • V6> 2 — £divvi — r/ 2 divw — /3divvi — 6> 2 divw, 

where J and / read as 

J={eQ + £)f - (eQ + £)(C/ 2 + v 2 ) ■ V(*7 2 + v 2 ) 

-(ePi + m)((W + w) • V(U 2 + v 2 ) + (Ui + vi) • V(W + w)) - V(M + Q0 2 ), 

J = 2// J D(W + w) : D(U 2 + v 2 ) + 2fiD(Ui + Vi) : D(W + w) 
+Adiv(W + w) • div(U 2 + v 2 ) + Adiv(C/i + vi) ■ div(W + w) 
-(eQ + 0(^2 + v 2 ) • V0 2 - (ePi + m )((W + w) • V0 2 + (Ui + vi) • V/3) 
+(eQ + £)0 2 divv 2 — (ePi + r/i)(/3divv 2 + 0idivw) + Qdivv 2 + Pidivw. 

Note that J and / can be bounded as follows. 

|o <(e||Q||i + ||e||i)(||f|| 2 + \\U 2 g + ||v 2 |||) + (c||fi|| 2 + |M| 2 )(||Wii + ||w||i) 

(3-11) 

2 + ||v 2 || 2 + ||C/i|| 2 + ||vi|| 2 ) + ||0|| 2 ||Q||i + ||/3||i||Pi|| 2 
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and 

||/||o <C{(||W||i + ||w||i)(||t7i|| 2 + \\U 2 \\ 2 + ||vi|| 2 + ||v 2 || 2 ) + (e||Q||i + U\\i)(\\U 2 \\ 2 + ||v 2 || 2 ) 
• \\O2W2 + + lhi||2)[(||W||i + ||w||i)|N| 2 + Ma + l|vi|| 2 ||/3||i] 

+ (e||Pi|| 2 + hi||2)(||^||i||v 2 || 2 + ||#i|| 2 ||w||i) + ||Q||i||v 2 ||2 + ||Pi||2||w||i} 

(3.12) 

On the one hand, we multiply (|3.9p i by W and make use of Poincare's inequality to deduce 
that 

^-9Mi)\\ VW f < C (\\w\\i + H|i)||c/ 2 || 3 ||w||i, 



2 

where C is a positive constant depending only on Q, and fi. Consequently, 



(3.14) 



|i<C(||W||i + ||w||i) (3.13) 

for some positive constant C depending only on fi, //, A, f , E and eo- 
By the classical estimates for the Stokes equations 

j —fj,AW + VQ = -{W + w) • VU 2 - {Ui + v x ) • VW, 
\ divW = 0, 

we obtain that 

||Wi 2 + ||VQ||o <C(]|T^ + w||i]|C/ 2 ]| 2 + + vi]| 2 ||VF||i: 

<C(||W||i + ||w||i), 

where the estimate (|3.13p has been used. 

On the other hand, if we multiply (|3.10p i . (|3.10p 2 and (|3.10p 3 by £, w and /3 in I? respectively, 
we find that 

H|Vw||2 + C||divw||2 + K||V/3||2 

= - j £divvi + r/ 2 divw + v 1 ■ V£ + w • Vr/ 2 + ediv(Pi(iy + w) + Q(U 2 + v 2 )) + W ■ Vr] 2 £dx 
+ J eJ - w • Vv 2 - vi • Vw - W • Vv 2 - V(0i£ + h: 



wdx 



+ 



eD — w • V9 2 — vi • V/3 — W ■ V6 2 — £divvi — r/ 2 divw — /3divvi — #idivw (3dx 



< C{||^|| 2 |hl|3 + [|fl|ill*llill»»lli + eUhiWPiHWWh + HwllO + ||Q||i(||C/ 2 || 2 + ||v a || a )] 
+ l|W / llih 2 ||3||el|i} + C{e\\ J\\l + ||w||f(||v 2 ||l + Hvilli) + ||W||?||v 2 ||I + 

Mil} + c{e\\i\\i + iiwiif n^iii + iiviinnfli! + \\w\\i\\e 2 \\i + neii?i|viiil 

+ IM||||w||? + ||^||?||vi||I + ||^i||i||w||?} + f(||w||g + Ml). (3.15) 

Also, from the Stokes equations 

j -fie Aw + V£ = e(eJ - w • Vv 2 - vi • Vw - W ■ Vv 2 - V(0i£ + f3r] 2 ) + (Vdivw) - V/3, 
1 divw = divw, 

we get the following estimate 

e||w|| 2 + ||V£|| <eC( ||divw|| 2 + e|| J|| + ||w||i(||v 2 || 2 + ||vi|| 2 ) + ||^||i||v 2 || 2 



(3.16) 

-HKMblKILi -UMIklMh) + <T'||-*||j. 
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Applying Poincare's inequality and substituting (|3.16|) into (|3.15p . employing the estimates 
(|3.1ip ~ ()3.14p and recalling the smallness of eo and E, we conclude that 

+ ||w||i + ||/3||i < C{\\W\\! + ||w||i + ||/3||i), 

where C is a positive constant depending only on A,f, E and eo- This completes the 
proof. □ 

Finally, having had Lemmas 13.11 and 13.21 we can apply the Tikhonov fixed point theorem to 
find a fixed point (U,v,9) = N(U,v,9) in the set K. Moreover, the pressure P G H 2 satisfies 

VP = f + g + fiAU - (U + v) • W, 

and (U + v, eP + n, 9) is a solution to (11.5p . Thus, we have shown the following proposition. 

Proposition 3.1. Let f, g € H 2 (Q,). Then, there exists an eo depending only on $7,/i,A,f and 
g, such that for all e G (0,eo), there is a solution (U,P,v,r/,9) £ (H A n Hq a ) x H 3 x (i7 3 n 
Fq 1 ) x # 2 x (# 3 n H%) of (J23]) and satisfying 

l|v|| 3 + + ||#|| 3 < 

where E is a small positive constant depending only on 0,yu,A, f and g. Moreover, (U + v, eP + 
n,#) is a solution of the system (|1.5p /or any e £ (0, eo). 

4 Incompressible limit 

Let e < eo and (t/ e , v e , 6* e ) G -ftT be the solution established in Proposition 13.11 We take 
v = v = v e , 9 = 9 = 9 e and n = rf in (j3.6|) to get that 

K|| 3 + ll^lls + Mh <Ci 6 (l + M!) 5 {e[||f ||i + (Mi + l) 2 ] + (E + E$) + <5}||n|| 2 

+ (eM? +E)\\v\\ 3 + (5M? +E)\\9\\ 3 

+ Ci 6 (l + Mi) 5 {e(Mi + l) 2 + e(l + Mi) [||f ||i + (Mi + l) 2 ] }. 
Thus, by taking eo and E so small that 

Ci 6 (l + Mi) 5 {e[||f||i + (Mi + l) 2 ] + (P + pi)} < 1, eM? + E < 1, 

we obtain 

l|v e || 3 + Wh + lh £ ||2 < Ci 6 (l + Mi) 5 {e(Mi + l) 2 + e(l + M x ) [||f ||i + (M 1 + l) 2 ] }, 
whence, 



l|v e || 3 + \\9 e h + \\v e h -> °> as e ^°- 

Furthermore, from (|2.2p i . i.e., 

divv e 



(4.1) 



and ()4.ip we get that as e — > 
divv 



-v e • Vrf - n e divv e - ediv(P e (U e + v e )) 



<||v e • Vn £ ||i + ||ndivv £ ||i + ||ediv(P e (£T + v e ))||i + \\U e ■ Vn e ||i -> 0. (4.2) 
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Due to (03]) and 



Vrf + V0 e 



eF e - v e • Vv e - FVrf - ?f W - U e ■ Vv £ + //Av e + (/z + A)Vdivv £ 



with F € = (eP e + rj e )f - (eP e + rf)([/ e + v £ ) • V(C/ e + v e ) - 9 e VP e - P £ V6 e , which comes from 
the transform of (|2.2P o. one deduces, recalling Poincare's inequality, that 



rf + 



0, 



as e 



0. 



(4.3) 



On the other hand, in view of Lemma [2 .21 we observe that (U e , P e ) is a uniform-in-e bounded 
sequence in (H 4 n Hq) x H 3 . Hence, there are a subsequence of (U €h ,P 6k ), still denoted by 
(U 6k , P 6k ) for simplicity, and (U, P) £ (H 4 n Hq) x H 3 , such that as e -> 0, 



and 



(U e ,P e ) 



(U, P), weakly in (i/ 4 n H^) x fl" 2 



(C7, P), strongly in (i/ 3 DH^xH 2 . 

Thus, if we take to the limit as e — > in (|2.1[) and (|2.2[) . we conclude that ([7, P) is a solution 
of the steady incompressible Naiver-Stokes equations (|1.7j) . 
In conclusion, we have that 



lim inf ||£/ e +v <: 

e-5-Oi7,PeL 



Uh + 



P e + 



P 



+ 



0. 



where L is the same as in Theorem 11.11 Thus, the proof of the low Mach number limit is 
completed. 
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